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Abstract 

The general theory of the branching processes is used for establishing the 
relation between the parameters k and n of the negative binomial distribution. 
This relation gives the possibility to describe the overall data on multiplicity 
distributions in pp(pp)-col\isions for energies up to 900 GeV and to make several 
interesting predictions for higher energies. This general approach is free from 
ambiguities associated with the extrapolation of the parameter k to unity. 



1. Introduction 

Theoretical description of the multiple production, based on so-called soft processes, 
today is beyond the limits of QCD, and the natural approach there is to look for empirical 
relations. 

* On leave from Yerevan Physics Institute, 375036, Yerevan, Armenia 
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The most popular in this field was the KNO scaling for multiplicity distributions which 
was fulfilled very well for hadronic collisions up to ISR energies and for e + e _ -annihilation. 
The evident violation of the KNO scaling at energies of the CERN collider [1] attracted 
much attention to the negative binomial distribution (NBD) which describes fairly well the 
overall features of the data on the multiplicity distribution (MD) of hadrons in different 
processes (pp(pp), e + e~ , up, AA . . .), in different ranges of rapidity and in a wide interval 
of energies [2,3]. It is especially relevant to the pp(pp) interaction. 

Taking into account the special role of NBD in describing the multiplicity distribution 
at high energy, it seems to be important to consider NBD on the basis of general as- 
sumptions about the character of the process of particle production in hadronic collisions 
without detailed specification of the dynamics. In reference [4] it was proposed as a basis 
for NBD to consider the multiple production as a random stationary branching process 
which is a rather general probabilistic model for the processes of the multiplication and 
transformation of the active particles. In this approach the transformation of each particle 
is independent of the history of the process and of the transformation of other particles, 
obeying the general probabilistic laws of Markov processes. The same refers to the fate of 
the generation of each particle. 

The branching processes may find their realization in terms of the quark-gluonic cas- 
cades, corresponding to the microscopic description of the nonequilibrial evolution of the 
partonic system, e.g., in the rapidity space [5,6]. It is important to stress that for us there 
is no need to know the details of the dynamical laws governing these cascades. 

It was realized that the system of produced hadrons may be considered as a result of 
the contribution from coherent and chaotic components (so called two-component model) 
and it is known that in pp(pp) collisions at high energy the chaotic component [7-10] 
dominates, which described by NBD, whereas in e + e~-annihilation the coherent (Poisson) 
part is essential. 
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The observation of dynamical chaos in the dynamics of nonabelian gauge fields (see 
e.g. [11]) raises the question about the role and origin of the chaotic component in hadronic 
collisions. 

There exists an interesting practical observation [7,9], that in distinction to e + e~- 
annihilation where the addition of the small (~ 10-20%) chaotic (noise) amplitude es- 
sentially changes the multiplicity distribution, in pp(pp) -collisions the addition of a small 
coherent component to the NBD does not change the shape of the distribution significantly. 

Taking into consideration the above mentioned arguments and remarks we here con- 
sider the NBD as adequate for the description of the multiplicity distribution in pp(pp) 
collisions at high energy. 1 The NBD 



r(n + i)r(fc) + 



p (k) = u^ + fcj {k) m 



has two parameters n and k. n is the average multiplicity. As for k, initially it was asso- 
ciated with the number of chaotically emiting cells. After the UA5 experiments [1,10,14] 
it is clear that such a meaning of k in general is not necessary, because SppS collider data 
yield the empirical relation 

^ = a + 61nv^ (a « -0.1, b « 0.06) (2) 

which is valid up to energy ^fs = 900 GeV, not showing a tendency for saturation. So, at 
such energies the KNO-like scaling continues to be violated, which is more clearly expressed 
in the observed strong rise of the moments 



(n 



with y/s [14]. 



1 In principle on the basis of the quantum optics it is possible to generalize NBD to 
take into account the chaotic as well as the coherent components [7,12]. 
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Of course, one cannot extrapolate the concrete form of the empirical relation (1) to 
higher energy since this would lead to contradiction. For instance, from (1) it would follow 
that the peak of the distribution would be at n = at very high energy when k = 1. This 
means that saturation of k must take place at ultrahigh energies at a value larger than 
unity (see also [15]). It indicates the necessity to establish the relation between k and ^Js 
(or, at least, between k and n) based on general theoretical considerations. We propose 
that such a basis could be a general theory of branching processes. As mentioned such an 
approach was developed in [4] where the idea of the stationarity of the branching process 
was used for establishing the relation between k and n. The result 

- = a + b\n- (a « 0.12, 6^0.08) (4) 
k k 

gives the unconfined though weaker rise of ^ with ^Js leading to the above difficulty 
associated with the extropolation of k to unity. Unfortunately, in deriving (4) the authors 
of [4] incorrectly used the conditions for stationarity of the branching process. In the 
paper [16] the condition for stationarity is used correctly though the authors missed the 
most interesting ansatz, in our opinion, of the relation between k and n. The resulting 
dependence of k on n, 

k = A {^j B (A ~ 11, B « -0.5) (5) 
again did not avoid the problem resulting from the extrapolation of k to unity. 

2. Relation between k and n 

Thus, we consider NBD as a result of the stationary branching process with one sort of 
multiplied particles (pions) and continuous evolution parameter t. The generating function 
F for such a process satisfies the reverse Kolmogorov differential equation [17]: 

f = /(*«>■ (6) 



For the generating function of the NBD 

F(M) = ]TPi fc V= [l + £(l -a)]"* (7) 

n 

f(F,t) equals 

-FlnFy + F(l-F 1/k )k — , (8) 
k m 

where m = n/k and k = dk/dt, etc. 

For a stationary branching process /(F, t) is factorized, /(F, t) = (p(F)ip(t). Evidently 
the condition 

k . m , nS 

— = const, k — (9) 
k m 

[4] which leads to the relation (4) does not give such a factorization. For F ~ 1 factorization 
takes place [16] and this approximation is also good for k close to unity. More adequate 
here is a parameter 

s = — (10) 

which is small at F ps 1 and not too small k. Expanding 1 — F x l k up to S 2 in (8), it is 
easy to obtain the solution of the resulting differential equation which is a necessary and 
sufficient condition for factorization: 

n — o 

where a and b are the integration constants which one must find from comparison with 
experimental data. Thus it is possible to state that NBD with relation (11) between its 
parameters k and n is the consequence of a stationary branching process. 

The function k(n) is very simple. At ab > k is decreasing from a to — oo and from 
+oo to a. But experimentally (at least for 10 < yfs < 900 GeV) k is decreasing with ^Js 
[1,14,10], so physically interesting is a case ab > 0. But at the same time the case a < 0, 
b < is also unphysical, because it corresponds to n < 0, or to k < 0. By the same reason, 
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if we do not want to have negative k, we must discard the lower branch of (5) with ab > 
corresponding to decreasing k in (a, — oo) interval. 2 

Thus this simple analysis has shown that in the region of high energy the relation 
between k and n is given by (5) with positive a, 6, and it is necessary to confine to the 
branch of the hyperbola (5) in the first quadrant. 

3. Consequences and predictions of the model 

The relation (5) between k and n, in spite of its simplicity, is rather rich in content. 
Let us stress once more that this relation must be only used for n > b p=j 7, i.e. for y/s > 14 
GeV; smaller energies should not be considered here. Our model ensures that k > a > 
(from the fit follows a = 3.06), implying that the limit k = 1 never is achieved. Thus there 
does not exist the difficulty associated with k = 1 at very high energy that is characteristic 
of some other ansatze [1,4,16]. 

From (5) it follows that asymptotically, when k goes to the saturation, KNO scaling is 
restored. The asymptotic distribution function at y/s » 14 GeV and n 3> k = a ps 3.06 
has a form of T-distribution: 

i;(z) = fiPW » kkzk ~[ e ~ kZ = UA9z 2m e- 3mz (z = -) (12) 

Our model gives rather clear predictions for C g -moments. In particular, Wroblewski's 
relation here takes place well at high energies: 



D 2 Vn 2 -n 2 . i (\ IV ( 1-92 \ 

~ ' (C 2 -l)* = ( T + -) w0.57 1-— . (13) 



n n \k n ) V n 



2 In this connection there is an interesting observation in [18] that the multiplicity data 
for small (y/s < lOGeV) energies is possible to describe well with negative k. From this 
point of view it may be said that two branches of the hyperbola (5) naturally divide the 
large energy region (n > b > 0) from small energies (n < b). From our fit (see below) 
b ps 7, so it means that we must not consider energies below yfs ps 14 GeV. 
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The second order correlation g^> = ^p-^ which is increasing slowly and asymptoti- 
cally equals 1.33 indicates not the presence of a coherent component as sometimes stated 
but just the fact that k is always larger than unity. All high order moments C q are rising 
and saturate asymptotically, as is easy to understand from the relation (q > 2) 

q ~ 2 /1\ 

c q = i + £ k (c 2 - ly-™- 1 (14) 

m=0 ^ ' 

where P^ (^) are polynomials of order m with positive coefficients (P = 1) and ^ = 
is increasing. Asymptotically we have: 



1 _ bl 

a an 



C ~ T{k + q) (15) 



i.e. (up to 1/n 2 ): 



1.27 

C 2 « 1.33 — 

n 

C 3 « 2.21 - — 
n 
21 3 

C 4 « 4.39 

n 

C 5 ^ 10.19-^^. (16) 

n 

The scaled peaks of the multiplicity distributions is moving to the left toward its 
asymptotic value: 



%eak(Vs) = — = 1 h-- = 0.67 + — . (17) 

n a an n 



Finally, before going into the comparison with experimental data, let us make one 
comment. By no means do we consider the limiting value of k = a as an indication that 
corresponds to asymptotic value of the number of clusters, fireballs, minijets etc. in the 
multiple production. 3 



3 Note that some experiments (see [20]) are indicating that at sufficiently high energies 
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Let us stress only that the often used value of k = 1 is meaningless. 4 In particular, in 
connection with this value of k in [16] it was made very strong and unusual statement that 
in the process of the multiple production the information entropy achieves its maximal 
value for k = 1 and as a result n achieves its maximal value and thus does not depend on 
the energy at all. This statement is derived from the fact that this entropy for NBD near 
k = 1 behaves as Inn + 1 — — lj (k — l) 2 . But the lower bound on k > a in our model 
shows that such a statement is a result of the unphysical interpolation of k to unity. 

4. Comparison with experiments 

To obtain numerical values of constants a and b in (5) we used the results of the fit of 
parameters n and k of NBD by experimental distributions of charged particles multiplicity 
in the range from y/s = 19.5 GeV to \fs = 900 GeV [10,14,19] (non-single diffractive 
events). The result of the fit of a and b in (5) on the basis of these data gives: 

a = 3.06 ±0.06, b = 6.95 ± 0.08. (18) 

In Fig. 1 is shown the function (5) obtained with these values for a and b. We did not 
consider points corresponding to low energy (see footnote 2). Fig. 2 shows the dependence 
of 1/k on n, which is seen to saturate. Fig. 3 gives the curves for C q (q = 2 — 5) for our 

(E L ~ 400 GeV) the number of the clusters produced in pp-interactions is 4.2 ± 1.7. If 
we continue such an interpretation of /c, then k~ x may be considered as the ratio of the 
probability for two particles to be emitted from one cluster to the probability of emission 
of these particles by two different clusters [6]. So the asymptotic "aggregation" degree is 
a" 1 « 0.33. 

4 The k = 1 in our model corresponds to the negative n (lower branch of (5) from a 
to — oo). Notice that k > a ~ 3 shows that our expansion parameter of 5 = k~ x InF is 
adequate and selfconsistent. 
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model (solid line) and compares them with experimental data [10,14,19] for non-single- 
diffractive component of pp(pp) reactions. It is seen that higher moments (q = 4, 5) have 
not yet achieved their asymptotic values (C q = 4.39 and C5 = 10.19 at yfs = 900 GeV. 

In Fig. 4-7 are shown the distributions ip(z, k) as a function of z = n/n for energies 
0.546, 1.8, 8, and 40 TeV, respectively. Solid lines show the asymptotic distribution (12). 
It is seen from these figures that scaling sets in at y^s 8 TeV. These curves show the 
systematic shift to the left of the peaks of ip(z, k) with increasing y/s. 

Finally, Fig. 8 shows the information entropy 



which is defined by the chaotic component only for k from (5). The figure also shows the 
"maximal" entropy iu max corresponding to k = 1 which is meaningless in our model. 

5. Conclusions 

In the present paper we have attempted to establish the relation between parameters 
of k and n of NBD on the basis of the general theory of random branching processes. 
This relation seems to be rather interesting, selfconsistent and has a predictive power. It 
removes some contradictions which occured in the use of NBD for description of multiplicity 
distributions for high energy hadronic collisions. 

On the whole the agreement of our model with the existing experimental data is good 
enough which, of course, is not surprising because of the coefficients a and b in (5) were 
derived from a fit to the experimental data for k and n. More important are the predictions 
for the behavior of C q (^-) and ip{z, k) at higher energies which may be checked at LHC 
and SSC: the restoration of the KNO scaling in the multi-TeV region, asymptotic constant 
values of C q , depending on q, "explanation" of the Wroblewski rule at high energy, and 
the asymptotic value of the peak of ip{z). 




(19) 
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It is interesting to compare qualitatively our model of general branching processes 
with the detailed models of quark gluonic branching processes. If one neglects the quark 
branching the resulting parton distribution looks very similar to NBD and their conclusions 
qualitatively coincide with ours (limit of the widening of the distribution shape, increase 
and final saturation of C q , etc. [21] .) The dominant role of gluonic branching in comparison 
with quark branching is the characteristic feature of the detailed study of corresponding 
processes from the point of view of dynamical chaos [22] , or from the approach based on the 
detailed consideration of branching of quarks and gluons at the formation of quark-gluon 
plasma [23,24]. 

There is at least one aspect which apparently necessitates the quark branching: the 
observed small oscillations in the high-multiplicity tail of P n -distribution at Tevatron en- 
ergy [25]. If we recall the very old prediction of such oscillations in the Regge-pole approach 
[26] which is connected with Pomeron cuts, then it seems reasonable that quarks may be 
responsible for these phenomena. (The "explanation" of these oscillations by the addition 
of two binomial distributions (five-parameter fit [25]) may also be the reflection of this 
two-Reggeon cut. 

Finally in connection with the meaning of parameter k of NBD and its asymptotic 
limit in our model (fc m in ~ 3) it would be very interesting to apply our model to the 
multiplicity distribution of hadrons in irp, as well in e + e~, up and ep collisions at high 
energies. 

In conclusion we thank to I. G. Aznaurian and N. L. Ter-Isaakian for discussions. One 
of us (S.G.M) is grateful to Berndt Miiller for interesting discussions and useful advice, 
and to C. Gong and C. Wang for help with the fit of the experimental data. S.G.M. thanks 
the Department of Physics, Duke University for hospitality. This work was supported in 
part by the Department of Energy (Grant DE-FG05-90ER40592). 
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Figure Captions: 

Fig. 1: The dependence of k on n from eq. (5) with coefficients a = 3.06 b = 6.95 obtained 
by fit. Two points shown on Fig. 1 and corresponding to low energies < 19.5 
GeV), are not taken into account (see footnote 2). 

Fig. 2: -i dependence on n from (5) (a = 3.06, b = 6.95). 

Fig. 3: The C^-moments (q = 2 — 5) as a function of n from (5) (solid lines) compared 
with experimental data on inelastic, non-single-diffractive component of pp(pp) 
reactions (see table 2 from [27] and [10,14]). 

Fig. 4: The dependence of nP n on z = Dashed line for = 546 GeV, solid line is 
an asymptotic distribution. 

Fig. 5: Same for = 1800 GeV. 

Fig. 6: Same for = 8 TeV. 

Fig. 7: Same for = 40 TeV. 

Fig. 8: The information entropy as a function of n. Solid line is for our model with 
k dependence of (5). Dashed-dotted line corresponds to the "maximal" entropy 
(k=l). 
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